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NuTeV is a neutrino-nucleon deep-inelastic scattering experiment at Fermilab. The detector consists of an iron-scintillator sampling 
calorimeter interspersed with drift chambers, followed by a muon toroidal spectrometer. We present determinations of response 
and resolution functions of the NuTeV calorimeter for electrons, hadrons, and muons over an energy range of 4.8 to 190 GeV. 
The absolute hadronic energy scale is determined to an accuracy of 0.43% . We compare our measurements to predictions from 
calorimeter theory and GEANT3 simulations. 



1 Introduction 

The increased intensity of the Fermilab Tevatron fixed- 
target program has made it possible to qualitatively 
improve neutrino deep-inelastic scattering experiments. 
Deep-inelastic neutrino scattering probes both the elec- 
troweak and strong forces in unique ways which are both 
competitive and complementary to other measurements 
at hadron and electron colliders. For these reasons, it is 
important to continue improving the precision of mea- 
surements with neutrino beams. NuTeV (Fermilab Ex- 
periment 815) is designed to exploit the intensity capa- 
bilities at Fermilab using a new neutrino beam, an up- 
graded neutrino detector, and a continuous test beam 
calibration system. 

The new neutrino Jpcam uses a sign-selected, 
quadrupole train (SSQT) 111 to produce a high-intensity, 
ultra-pure beam of cither neutrinos or antincutrinos. For 
neutrino detection, the experiment uses an upgraded ver- 
sion of the CCFR detectorB, shown in Figure I, with new 
scintillation oil and photomultiplier tubes combined with 
refurbished drift chambers. In the detector, neutrino in- 
teractions produce a hadronic shower from the outgoing 
struck quark whose energy is measured in the target- 
calorimeter and, for charged current events, an associated 
outgoing muon whose angle (momentum) is measured in 
the target-calorimeter (downstream muon spectrometer). 

The NuTeV data run took place during the Fermi- 




Figure 1: The NuTeV neutrino detector showing the target 
calorimeter followed by the downstream muon spectrometer. 

lab 1996-1997 fixed target run. The experiment recorded 
over three million neutrino and antineutrino interactions. 
Two of the physics goals of NuTeV are a precise measure- 
ment of the weak mixing angle and measurement of struc- 
ture functions and the strong coupling constant from 
QCD scaling violations. Both of these results are depen- 
dent upon a detailed understanding of the response of 
the target-calorimeter. A previous experiment using this 
calorimeter, CCFR, determined the calorimeter energy 
scale to an uncertainty of approximately 1%. CCFR mea- 
suresB the weak mixing angle to be sin^ 9w = 0.2236 ± 
0.0019(stat) ± 0.0019(syst) ± 0.0030(model). NuTeV 
aims for a total precision of better than 0.002 on sin^ 9w, 
primarily by changing the measurement technique to re- 
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duce model uncertainties. However, in CCFR the ex- 
perimental systematic uncertainty due to calorimeter re- 
sponse is ±0.0011 and the NuTeV technique is consider- 
ably more sensitive to energy calibrations. The reduced 
theoretical uncertainties make an improved calibration 
essential for the success of this measurement. 

In the case of the measurement of the strong cou- 
pling constant, the systematic uncertainty on the QCD 
scale parameter, A-p^, frorrLcalibration effects is at the 
50-100 MeV level previousljQ, which is the largest single 
experimental source of uncertainty in the measurement. 
In NuTeV, this uncertainty would be reduced by a factor 
of three by an absolute calibration of 0.3% uncertainty. 

For the measurements listed above, the absolute re- 
sponse and resolution of the hadronic shower energy mea- 
surement are crucially important. For this reason, pre- 
cision detector calibration and response determination 
is a key component of the NuTeV program. This is 
done using several data sets: the actual neutrino events, 
neutrino-induced muons from upstream shielding, and 
calibration beam data. Throughout the data run, the 
calibration beam operates continuously and provides mo- 
mentum tagged electrons, muons, and hadrons with en- 
ergies between 4.8 and 190 GeV. A precision spectrome- 
ter provides an event-by-event momentum determination 
with resolution better than 0.3% and a combination of a 
Cerenkov counter and a TRD are used to determine the 
particle type for each event. 

This article describes the various techniques and 
studies that lead to the precision calibration of this 
calorimeter. First, the detector and electronics calibra- 
tions using neutrino-induced events are described, fol- 
lowed by the test of these techniques and resolution stud- 
ies using the calibration beam data. This article also 
demonstrates which aspects of the detector response are 
accurately modeled by GEANT and other software sim- 
ulation packages. 

2 The NuTeV Calorimeter 

The NuTeV calorimeter consists of 168 plates of steel 
measuring 3 m (H) x 3 m (W) x 5.1 cm (L), interspersed 
with 84 scintillation counters of dimension 3 m (H) x 

3 m (W) X 2.5 cm (L) and 42 drift chambers. There are 
two plates of steel between every two consecutive scin- 
tillation counters, and one drift chamber between every 
other set of counters. One unit counter consists of a 
scintillation counter and two steel plates surrounding the 
scintillator. Therefore one unit calorimeter layer consists 
of two counters and a drift chamber. This configuration 
leads to a detector with 10.35 cm of steel between coun- 
ters and 20.7 cm of steel between drift chambers. The 
geometry of one unit of the calorimeter is shown in Fig- 
ure and this unit is repeated 42 times to make up the 
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Figure 2: Geometry of one unit of the calorimeter. This unit is 
repeated 42 times to make up the entire calorimeter. One unit 
of the calorimeter consists of a scintillation counter sandwiched 
between two steel plates. 



entire calorimeter. Table |l| summarizes the materials and 
their longitudinal sizes in units of cm, radiation length, 
and interaction length, for one unit of the calorimeter's 
longitudinal layer. 

Figure ^ shows a schematic diagram of a NuTeV 
scintillation counter. The scintillation counters are In- 
cite boxes filled with Bicron 517L scintillator oil. The 
counters have 3 mm vertical Incite ribs spaced by 2.5- 
5.1 cm, depending on the lateral position of the ribs, and 
designed for structural support. Since these ribs do not 
scintillate, the counters are staggered so that the ribs 
are not aligned on the transverse plane throughout the 
length of the calorimeter. 

Each counter is surrounded by 8 wavelength-shifter 
bars, doped with green BBQ fluor, and is read out in four 
corners by photomultiplier tubes (PMT's), mounted one 
on each corner. The PMT's are 10-stage Hamamatsu 
R2154 phototubes with a green-extended photocathode, 
with gains set to about 10^. 

There is an air joint between the wavelength-shifter 
bars. The joints between the wavelength-shifter bars and 
the phototubes have 3 mm thick clear silicon jelly cookies 
for better optical and mechanical connections, as well as 



Table 1: Composition in interaction length and radiation length of 
one unit of the NuTeV calorimeter. This unit is repeated 42 times 
to make up the entire calorimeter. 



Component 




Length 






cm 


^0 


A/ 


4 Steel Plates 


20.7 


11.75 


1.24 


2 Scint. Counters 


13.0 


0.51 


0.16 


1 Drift Chamber 


3.7 


0.17 


0.03 


Total 


37.4 


12.43 


1.43 
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for PMT window protection. The cookies are made of 
Dow Corning Sylgard(R) 182 silicon elastomer and 182 
curing agent. Given this particular geometry and readout 
scheme, NuTeV observes muon signal distributions which 
are consistent with, on average, 30 photoelectrons for 
muons traversing the center of a counter. For muons 
closer to the edge of the counter, where light collection 
is more efficient, the number is higher. 



PMT 4 PMT 1 




PMT 3 PMT 2 

Figure 3: A schematic drawing of a NuTeV scintillation counter. 



3 Calorimeter Readout Electronics 

The design of the readout electronics is dictated by 
the requirement to accommodate a very large dynamic 
range of the signal using 11-bit analog-to-digital con- 
verters (ADC's). Minimum ionizing particle energy loss 
(MIP) in the calorimeter is approximately 0.15 GeV per 
unit counter while neutrino interaction induced hadronic 
showers could deposit up to 100 GeV into a single unit 
counter. Note that the actual energy deposit of a mini- 
mum ionizing particle in a single scintillation counter is 
approximately 4 MeV. 

Figure ^ shows a schematic diagram of the readout 
electronics system. The readout electronic channels con- 
sist of the following three separate gains, to measure en- 
ergies in a wide dynamic range: 

• HIGH is the signal formed by a sum of signals from 
each of the four PMT's of a given counter (fan-in 
ES-7138) El and x 10 amplification of the summed 
signal by the linear amplifier LeCroy 612A. The 
amplified signal is then digitized by a LeCroy 4300Q 
Fast Encoding and Readout ADC (FERA). 



Counter 1 readout 




Figure 4: A schematic diagram of the NuTeV calorimeter readout 
electronics of a counter. 



• LOW is the signal from each of the 4 PMT's di- 
rectly digitized by LeCroy 4300 FERA. 

• SUPERLOW signals are the digitized sums of 8 
PMT signals which come from 8 different counters, 
each separated by 10 counters. Each PMTjSignal 
is attenuated by 1/10 (fan-in LeCroy 127FlIiI). 

A typical minimum ionizing particle signal produces 
80 ADC counts in the HIGH channel, 2 ADC counts 
in each LOW channel, and ~ ADC counts in each 
SUPERLOW channel. The LOW and SUPERLOW 
channels are calibrated with respect to the HIGH chan- 
nels. 

Hadronic showers from neutrino interactions deposit 
up to 600 GeV in the calorimeter (with maximum en- 
ergy deposition in a single counter of about 100 GeV). 
A typical hadronic shower in the calorimeter saturates 
the HIGH channels and leaves a signal of a few hun- 
dred ADC counts in the LOW channels. Thus the LOW 
channels are used to measure the shower energy. In a 
very small fraction of the time, one of the four LOW 
channels of the counter is saturated when the transverse 
position of a neutrino interaction is close to one of the 
PMT's. In these cases, the attenuated SUPERLOW 
channel is used. 

Calibration of the NuTeV calorimeter begins with 
the calibration of the readout electronics and tying the 
HIGH, LOW, and SUPERLOW channels into a hnear 
model. 
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4 Pedestal Subtractions 

The best way to determine the pedestal values of the 
ADC channels well is measuring them under exactly the 
same conditions as the neutrino data. To achieve this 
we use two different methods - one using a specially de- 
signed random trigger and the other using "quiet regions" 
of the calorimeter during real neutrino events. A random 
trigger is activated throughout the run in all gate types 
to measure pedestals under the same condition as the 
trigger of interest. The rate of the pedestal trigger is 
prescaled to provide the necessary number of pedestal 
events - typically 10 events per accelerator cycle - with- 
out overloading the bandwidth. 

The second method uses events where the trigger 
(T2) is designed for neutral current interactions and re- 
quires significant energy deposition in consecutive coun- 
ters in the calorimeter. For each T2 event, an offline 
analysis program finds "quiet regions" in the calorime- 
ter, using the following algorithm: 

1. Count the number of counters with pulse heights 
more than 1/4 of a MIP (S-bit ON) and reject the 
event if this number is larger than 10. 

2. Select the counters with their own S-bits and the 
S-bits of their 4 closest neighboring counters OFF. 

3. Check that all three readout ADC channels of the 
selected counters described in Section ^ have en- 
ergies less than 0.3 MIP. The cut value is chosen 
to be much less than 1 MIP but much larger than 
the pedestal; for example in a HIGH ADC chan- 
nel, one MIP is ~70 ADC counts and a typical 
pedestal width is ~3 counts, thus the cutoff value 
of 20 counts is ~7 standard deviations from zero. 

4. Use the readouts of each ADC channel of the se- 
lected counters as the pedestal values. 

The offline analysis procedures use these pedestal 
events to keep a running average for each electronics 
channel, using both these methods. The two proce- 
dures for measuring pedestals agree to within .015 ADC 
counts in the LOWs, and .02 ADC counts in the HIGHs. 
This pedestal uncertainty would contribute a constant 
term of 32 MeV to the hadron energy resolution if all 
pedestal differences were correlated, and a 3.6 MeV width 
if these pedestal differences were uncorrelated. The con- 
stant term in the hadron energy resolution is consistent 
with zero with an error of 110 MeV (see Figure The 
"quiet region" method is used for neutrino data pedestal 
subtraction. 

Because of the differences in the upstream magnet 
currents and detector environment, the pedestals during 
the calibration beam gate and the neutrino beam gate 



are not necessarily equal. In fact, some channels differ 
by as much as 0.3 ADC counts in the LOWs. For the 
analyses of the calibration beam data, we use the ran- 
dom trigger method to measure the pedestals because 
the upstream part of the calorimeter, where the calibra- 
tion hadron beam enters, always has energy deposited 
in every event, precluding the "quiet region" method. 
The neutrino data pedestal comparisons ensure that this 
treatment is completely accurate to the few MeV level. 

5 Electronics Cross Calibration 

Relative calibration of the different channels of electron- 
ics is needed because the minimum ionizing particle sig- 
nal is measured with the HIGH channels, while the neu- 
trino interaction signal is measured with the LOW and 
SUPERLOW channels. Following the assumption of 
linearity of all the components of the readout electronics, 
we assume that the HIGH channel is the linear combi- 
nation of the 4 LOWs of the same counter: 

4 

HIGH{i) ^^Rf{i) X LOWj{i), (1) 
i=i 

where i is the counter index, Rj\i) is the relative cali- 
bration constant between the LOW signal of the PMT j 
of the counter i and the HIGH signal of the counter i. 

The SUPERLOW is the linear combination of the 
8 LOWs, as follows: 

8 

SUPERLOW{i) = -Rj' (*) X I^OWj (k), (2) 
fc=i 

where i is the SUPERLOW channel number, k is the 
counter index, j is the fixed PMT index, and Rfii) is the 
relative calibration constant between the SUPERLOW 
channel i and the LOW channel of the PMT j of the 
counter k. The set of calibration coefficients R^^{k) and 
Rj^k) is calculated for every data-taking run using the 
least squares method. 

6 Counter Gain Position and Time Dependence 

As do all high energy neutrino experiments, NuTeV 
has an ideal calibration source to track counter gains: 
muons coming from upstream neutrino interactions in 
the shielding. This means that not only are the muons 
completely correlated in time with the actual neutrino 
beam, but also illuminate the detector in a similar fash- 
ion. This section describes how the sample of muons 
traversing the entire length of the detector during the 
neutrino beam is used to monitor the position and time 
dependencies of the individual counter gains. 
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Figure 5: Typical energy deposition of muons traversing one scin- 
tillator counter in units of ADC counts. 



Figure ^ shows a typical energy deposition profile for 
muons traversing a particular counter. There are, on av- 
erage, 30 photoelectrons per MIP per counter. Events 
with very low pulse heights are from particles that go 
through the ribs of a counter while the events with large 
pulse height are from muon bremsstrahlung and e^e~ 
pair production. Since pair production increases with in- 
creasing muon energy, this energy deposition pattern can 
be used as an event-by-event muon momentum measure- 
ment, as described in Reference 7. 

In order to characterize the distribution shown in 
Figure |^ in a way that is stable with respect to cuts, we 
use the truncated mean procedure^. The truncated mean 
is determined by calculating the mean of the distribution 
using all events, then taking the mean again but only 
including the events between 0.2 and 2 times the previous 
mean. This procedure is iterated several times until the 
difference between the mean of two consecutive iterations 
is less than 0.1% of the previous mean. Corrections are 
made event- by- event for the muon's momentum as well 
as the angle with respect to the direction perpendicular 
to the counter. This procedure provides a "mean" that 
is insensitive to the width and tails of the pulse height 
distribution. The truncated mean for 77 GeV muons is 
defined as 1 MIP. 




Figure 6: Average counter response to muons traversing as a 
function of position in the counter. The coordinates are normalized 
to the half width of the counter, 1.5 m, on both axes. 

6.1 Position Dependence 

We calculate the average position dependence of the trun- 
cated means (which will be referred to as the muon re- 
sponse) for each counter, by averaging over the entire 
neutrino run. Figure ^ shows that response as a function 
of position for a typical counter. The light collection is 
largest at the corners of the counters, near the four photo- 
tubes, as expected. The technique used to track the time 
dependence of the gains alters the position dependence 
of each counter as a function of time. This is because 
the gain for each phototube is determined independently 
as a function of time and is reflected into the position 
dependence. 

6.2 Time Dependence 

The gain of a single counter at a particular moment in 
time during the run depends on the gains of the four 
phototubes as well as the gain of the scintillator oil itself. 
To determine the time dependence of the counter gains, 
we calculate a fractional phototube map, which is defined 
as the fraction of light reaching a given phototube as a 
function of position within an independent counter. 

These fractional phototube maps are measured using 
high energy neutrino interactions where the pulse height 
is high enough to be seen in an individual LOW channel. 
Figure |^ shows two sample phototube maps. Notice that, 
as expected, the fractional maps are strongly peaked near 
the phototubes themselves and drop off sharply where 
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two wavelength-shifter bars meet in the center of the 
counter. 

We then fit the muon response over a short period of 
time to a function with four parameters, where the pa- 
rameters correspond to the gains of the four phototubes. 
So, if M^{x, y, i) is the run-averaged muon response map 
for counter i as a function of x and y, and Fj(x,y,i) is 
the fraction of hght reaching phototube j of counter i, 
then the run-averaged muon response for phototube j of 
counter i, Pj{x,y,i), is simply 

Pj {x, y, i) = M\x, y, i) x {x, y, i). (3) 

The time-dependent function is then 

M{x,y,i,t) ^J29jiht)Pjix,y,i), (4) 

3 

where gj {i, t) is the relative gain of phototube j of counter 
i at time t compared to the average gain over the entire 
run. 

Gains vary by as much as 10% as a function of time 
and vary, on average, — 0.16%/°C (compared to the value 
-0.n%/°C specified by the vendorEfl), as a function of 
temperature as shown in Figure ^. The temperature de- 
pendence varies from counter to counter by ±0.04%/°C 



7 Hadron Gain Balance 

The technique described in the previous section deter- 
mines the gain for each scintillator counter relative to all 
the others, using muons that traverse the entire detector. 
To set the absolute hadron energy scale of the detector we 
measure the response of the calorimeter to a monochro- 
matic beam of hadrons incident on the most upstream 



16 17 18 19 20 21 22 23 24 25 26 
Upstream Cart Temperature (°C) 

Figure 8: Sum of phototube gain coefficients, normalized to tlie 
average gain over the entire neutrino run period, as a function of 
temperature of the most upstream region of the calorimeter. 

part of the calorimeter (see Section |2| for details) . Since 
this hadron beam deposits all of its energy in the most 
upstream 10-12 counters, the hadron response of only 
those counters is measured. The hadron response mea- 
sured this way would be usable for the entire calorimeter 
if the hadron response were completely correlated to the 
muon response. However, geometric non-uniformities in 
the calorimeter give rise to relative differences between 
the hadron and muon response, breaking the correlation. 

The NuTeV calorimeter measures the hadron energy 
by sampling the shower every 10 cm of steel. The en- 
ergy deposited by a hadronic shower in the scintillation 
counters is only a small fraction of the total energy de- 
posited in the detector. Therefore, variations in the pas- 
sive material surrounding each counter affect the average 
hadron signal sampled in that counter. In contrast, the 
muon signal is only dependent on the variations in the 
active material. Since the relative gain of a counter for 
hadrons may not be completely correlated with the rela- 
tive gain for muons, setting the hadron energy scale for 
the first 10-12 counters is not sufficient to set the scale 
for the entire detector. In this section we describe the 
technique used to measure the hadron/muon gain ratio 
for each counter, using neutrino interactions that occur 
throughout the entire calorimeter. 

Apart from the low interaction rates, neutrinos are 
a perfect relative hadron calibration source for the en- 
tire calorimeter. First of all, if the detector is far enough 
away from the neutrino production target, the energy dis- 
tribution of neutrinos interacting in the most upstream 
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counter of the calorimeter is the same as that of the neu- 
trinos interacting in the last counter. In a charged cur- 
rent Vfjt or 77^ interaction, which is the majority of neu- 
trino events that NuTeV sees, both a hadron shower and 
a muon are produced and deposit energy in the calorime- 
ter. As stated previously, the hadron shower deposits 
most of its energy in the first few counters after the 
neutrino interaction, while the muon deposits a small 
amount of energy in each counter over many counters, 
depending on its angle and energy. 

The average measured energy in the calorimeter from 
neutrino interactions should not depend on where the 
neutrino interaction occurs, assuming that one always 
measures the energy by summing over the same number 
of counters from the event vertex. If one sums over the 
first 10 counters after the event vertex, then the muon 
contribution to a 70GeV shower is about 3%. The ad- 
ditional muon energy deposited in the hadronic shower 
region would reduce the measured effect, but the amount 
by which the muon's presence changes the measurement 
is negligible compared to the statistical uncertainty on 
the gains. 

To determine the hadron/muon gain ratio for each 
counter, we need a sample of clearly identified neutrino 
interactions in the calorimeter. Events are selected by 
requiring a final state muon, visible by a minimum en- 
ergy deposition extending over at least 20 counters (2 m 
of steel equivalent) after the event vertex. To remove 
cosmic ray backgrounds, events are also required to have 
a reconstructed hadron energy greater than 20 GeV. The 
hadron energy is determined by summing over the ener- 
gies of the 10 counters following the event vertex. Since 
this energy cut ultimately depends on the relative gains 
obtained from the technique, the procedure must be it- 
erated. The average hadron energy of the events passing 
all cuts is about 70 GeV. 

Events are also required to occur at least 4 coun- 
ters from the upstream end of the calorimeter, 20 coun- 
ters before the downstream end of the calorimeter, and 
within 1.27 m of the center of the detector. These fiducial 
cuts ensure that the event is not induced by a charged 
particle entering from the side or front of the detector, 
and that the hadronic shower in the event is fully con- 
tained within the calorimeter. Because of these fiducial 
volume cuts, we are unable to use this technique to de- 
termine the relative gains of the 15 most downstream or 
4 most upstream counters. However, the most upstream 
counter hadron/muon gain ratios are determined using 
a similar technique, described here, with the calibration 
beam hadron data. The first four counters' gains are set 
by comparing the calibration beam hadron response of 
showers starting in the most upstream set of four coun- 
ters with those in the next set of four counters that are 
immediate downstream of the first set and whose gains 
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Figure 9: Average longitudinal hadronic shower profiles of neutrino 
events in two different shower energy ranges. It can be seen from 
the plots that large fraction of shower energy is deposited in two 
to three consecutive counters. 



have been determined from the neutrino data. 

The fitting procedure constrains the hadron energy 
of neutrino interactions to be constant by varying the rel- 
ative gains of the counters. Let the visible hadron energy 
of a neutrino event that starts in counter i be denoted as 
EHADi. In a given event the individual counter energies 
using the muon-derived gains are denoted by E{j), and 
the hadron/muon gain ratio for each counter is denoted 



by hj. In this notation. 



(5) 



EH AD, = i:]t\h,E{]) 

where the sum over 10 counters (1 m of steel equivalent) 
is expected to include more than 95% of the hadronic 
shower. The average hadron energy over all neutrino 
events that occur in counter i is then 



< EH AD, > 



^3 



i+9 



h,E{j) > 
hj < E{j) > 



Thus in theory the average hadron energy of showers that 
start in counter i depends not only on counter i's gain 
but also on the gains of the nine subsequent counters. In 
practice, however, hadron showers deposit a large frac- 
tion of their energy in only two or three consecutive coun- 
ters around the shower maximum, as shown in Figure ^ 

One first computes the average hadron energy 
{EAVE) over the entire fiducial volume by setting all the 
initial gains hj to unity, making the cuts described above, 
and calculating the average EHADi over all the events 
that pass the cuts. Then, one can fit for the hadron gains 
by minimizing a x^, where the is defined as follows: 

EAVE)^ 



X 



V80 



(£'iJAA(g»,g»-i,g»-2- 
errJ 



(6) 



and the ERRif is defined as: 



ERRf ^ {< EHAD^ > - < EHADi>^)/N, (7) 
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where Ni is the number of neutrino events, passing all 
cuts, that start at counter i, and the error is calculated 
assuming the gains hj are all set to their initial values 
(which for the first pass is simply unity). To fit for the 
gains, in theory, one simply has to minimize the de- 
fined in Eq. ^ determine the gains, and then iterate, 
making the energy cut and re-computing the new average 
and errors using the gains from the previous iteration. 

In fact the defined in Eq. ^ is unstable, due to 
the fact that a hadron deposits most of its energy in 
two or three consecutive counters at the shower maxi- 
mum and this causes a strong correlation between the 
two counters next to each other. The x^j when com- 
puted this way, is low not only for uniform gains very 
close to 1 but also for gains which are staggered by a cer- 
tain amount, where the even counters are all high and the 
odd counters are all low (or vice versa). This variation 
in hadron gains is larger and more regular than would 
be expected from detector non-uniformities in thickness 
and composition of material. This artificial hadron gain 
variation is avoided by separately fitting the gains using 
events whose showers start in every other counter (for 
example, even- numbered counters), then using the com- 
plementary set of events (for example, showers starting 
in odd-numbered counters) and refitting. The resulting 
gains for all counters are consistent between the two fits, 
have smaller errors than when all showers are included 
at once, and are much closer to unity. 

By averaging the two fit results and iterating, the 
gains are stable to better than 0.2% after 3 iterations. 
The statistical uncertainty on each relative gain is about 
0.9%, and is larger near the downstream edge of the de- 
tector where there are only events starting upstream of 
those counters. Figure ^ shows the gains obtained af- 
ter four iterations using the technique described in this 
section. The gains have an RMS of 2.3%, and are consis- 
tent with geometrical non-uniformities in the calorimeter 
(water bag thicknesses, steel plate thicknesses, etc.). 

The relative hadron gains allow a determination of 
the absolute hadron energy scale of the entire calorimeter 
by measuring the response of the upstream most counters 
to a monochromatic beam of hadrons. The total statisti- 
cal uncertainty on the hadron gains in the overlap region 
where there is both calibration beam and neutrino data 
is equal to 0.4% and dominates the overall uncertainty 
in the hadron energy scale. Each counter's hadron/muon 
gain ratio uncertainty (0.9%) is uncorrelated between the 
counters. The contributions of this uncertainty to the 
calibration beam energy measurement is reflected in the 
uncertainty in hadron response measurement and is neg- 
ligably small due to statistically random longitudinal de- 
velopment of hadron showers (see Section ^2|) . 

The relative gains obtained using the technique de- 
scribed in this section are used for the energy reconstruc- 
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Figure 10: Relative hadron/muon counter gains which arise from 
detector non-uniformities unrelated to scintillator thickness. 



tion in both the hadron and electron response measure- 
ments. 

8 The NuTeV Calibration Beam 

NuTeV is designed to include a simultaneous calibration 
beam separated from the neutrino beam by f .4 seconds, 
yet running within the same one- minute accelerator cycle 
(see Figure |T^). The calibration beam is used to set 
the absolute energy scale of the experiment, and also 
to measure the response of the calorimeter to hadrons, 
electrons, and muons, in order to properly simulate the 
detector. Finally, the calibration beam is instrumental 
in monitoring the time dependence as measured by the 
muon map technique described in Section |[ 

The calibration beam period within a cycle is f 8 sec- 
onds, and the typical beam incident angle to the center of 
the NuTeV calorimeter is 43 mrad in the horizontal direc- 
tion (0 mrad in vertical) with respect to the centerline of 
the calorimeter. The calibration beamline can transport 
particles of energies from 4.8 GeV to f 90 GeV, and de- 
pending on the beamline apparatus and magnet settings, 
can produce high purity beams of electrons, hadrons, or 
muons for energies above 30 GeV. 

The beamline is instrumented as a low mass spec- 
trometer with a long lever arm. The distance between 
the most upstream chambers in the spectrometer and 
the momentum-analyzing magnets is 83.3 m, and the 
distance between the most downstream chamber and the 
magnets is 69.2 m. This separation allows a modest align- 
ment uncertainty of 1 mm to translate into only a 0.1% 
uncertainty in the absolute momentum scale. The event- 
by-event resolution of the spectrometer, dominated by 
multiple scattering, is better than 0.3% for most energies. 
The beamline instrumentation is augmented for some of 
the run with a removable Gerenkov detector and a TRD 
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Figure 11: Accelerator time structure. Note that the interval be- 
tween the last neutrino ping and the slow spill calibration beam is 
only 1.4 seconds, allowing essentially an in situ calibration. 



array used to measure the beam particle composition. 

Over the course of the experiment, standard runs 
were taken at least once a week (50 and 100 GeV 
hadrons) and hadron energy scans between 4.8 GeV and 
190 GeV were taken once a month. Overall, NuTeV ac- 
cumulated a total of 17 million test beam triggers. 

8.1 Beam Time Structure 

The accelerator time structure during the 1996-97 Fer- 
milab fixed target run is depicted in Figure The 
accelerator complex cycles every 60.1 sec. The neutrino 
beam is delivered in five fast resonance extraction pulses 
("pings") of 5 msec width. The pings are separated by 
0.5 sec. The NuTeV slow spill calibration beam begins 

1.4 sec after the last ping and has a duration of 18 sec 
with uniformly distributed beam intensity. This calibra- 
tion beam is delivered in a beamline that is completely 
independent from the fast spill. This time structure 
provides continuous calibration data, taken concurrently 
with the neutrino beam, and allows for an in situ cali- 
bration of the detector. 

8.2 Beam Selection Scheme 

The NuTeV calibration program involves electrons, 
hadrons, and muons of momentum ranging from 4.8 GeV 
to 190 GeV. In order to select the desired type of parti- 
cles for a specific program with high purity, the beamline 
is designed as shown in Figure |l^. 

The target (NT8TGT) in the calibration beam is a 

7.5 cm (W) X 7.5 cm (H) x 30.3 cm (L) aluminum block. 
Protons of 800 GeV momentum strike the target with an 
integrated intensity between 4x10^^ and 8 x 10^^ protons 
on target throughout the 18 second long slow spill. The 
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Figure 12: NuTeV calibration beamline schematics. The "Ferris 
wheel" (NTACON) with four different thickness converter mate- 
rial is used to select pure hadrons or electrons. The 6 m long Be 
filter(NTBBE) is used to select pure muons. The numbers on the 
left-hand-side of each component indicate the relative distance of 
the component to the primary target (NT8TGT) in meters. Some 
beam position and intensity monitoring devices are not drawn in 
this figure because they are irrelevant for this paper. 



secondaries are then focused by a set of quadrupole mag- 
nets (NT9Q1 and NT9Q2) to the enclosure NTA and 
coUimated by a horizontal collimator (NT9GH) whose 
opening is adjusted depending on particle type and in- 
tensity. The polarity of the beamline is set to direct 
negatively charged particles to reduce intensity. 

The horizontal collimator (NT9GH) is then followed 
by a string of dipole magnets (NT9W-1, 2, 3, 4, and 5), 
fed by one power supply, for initial momentum selection, 
bending the beam in the horizontal plane. The verti- 
cal collimator (NT9CV) following the first set of dipoles 
(NT9W's) is used to further cut down the intensity. 

The "Ferris wheel" (NTACON), located immediately 
downstream of the set of collimators and the initial mo- 
mentum selection dipoles, has four mounts, each plac- 
ing a different thickness of material into the beam at 
a time. The thicknesses correspond to an empty hole, 
0.2Xo, 6Xo, and l2Xo. The empty hole is used for the 
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Figure 13: A schematic view of the NuTeV cahbration beam spec- 
trometer system. The large distances between the chamber stations 
allow accurate absolute momentum determination. 

muon mode, while the 0.2Xq piece is used for the electron 
niode. In the electron mode, the magnets downstream of 
the "Ferris wheel" are tuned for 20% lower momentum 
particles to compensate electron energy loss in the ma- 
terial. The two higher radiation length materials are put 
in the path of the beam for hadron modes to eliminates 
electrons from the beam. The thicker material is used 
for higher energy beams, the thinner for lower. 

The pinhole collimator (NTAPIN) following the 
"Ferris wheel" is put in the beam only for higher energy 
hadron modes {E > 30 GeV) to further cut down in- 
tensity and increase radiation safety. Typically, the size 
of the hole in the pinhole collimator is 5 mm x 5 mm 
and the momentum bite set by this collimator opening is 
approximately 0.2 GeV. 

The 6 m long beryllium filter (NTBBE) is only used 
in muon modes to filter out hadrons and electrons in the 
beam. The energy loss of muons in the filter is approxi- 
mately 6% while the survival probabilities of hadrons and 
electrons through the filter are less than 3.4 x 10^^ and 
3.5 X 10^^, respectively. The filter is then followed by 
additional two sets of dipoles (NTBWl-1, 2, 3, NTBW2- 
1, and 2) for further refinement of momentum selection. 
This combination of three large dipole strings throughout 
the long stretch of the beamline removes virtually all pos- 
sible contamination of unwanted particles and momenta. 

A final precision spectrometer is used to determine 
the beam momentum on an event-by-event basis. The 
spectrometer begins with two small area drift chambers 
(SDWCl and SDWC2 in Figure H) positioned at the 
downstream end of the last dipole string in the same 
beam enclosure (NTB). Figure [l^ shows a schematic view 
of the NuTeV calibration beam spectrometer system. 
The particle ID system, which consists of a Cerenkov 
counter followed by an array of TRD's, is located just 
upstream of the second set of chambers which, in turn, 
were positioned immediately upstream of the spectrome- 
ter dipole magnet string. When particle identification is 
not needed, these detectors are rolled out of the beam- 
line and are replaced by a vacuum pipe to reduce multiple 




Figure 14: PMT signals of the Cerenkov counter with a 160 torr 
nitrogen gas. A clean particle separation between 77(01), K(C2), 
and p(ped) at 50 GeV is apparent. 

scattering. 

The last dipole in the spectrometer magnet string 
can be rotated. This dipole is an integral part of the spec- 
trometer for beams with energies greater than or equal 
to 120 GeV and is also used to direct the beam to various 
positions on the detector surface for position dependent 
response measurements. 

8.3 Particle Identification and Beam Purity 

The Cerenkov counter provides particle identification for 
pions, kaons, anti-protons, and electrons, depending on 
the type of gas and the threshold pressure for each type 
of particles. The Cerenkov counter is equipped with two 
PMT's, Ci and C2, placed to face opposite directions, 
and is designed to act as a differential Cerenkov counter. 
The PMT Ci, accepts low angle (< 4.5 mrad) Cerenkov 
light from heavy particles, while the PMT C2, accepts 
large angle light from lighter particles of the same mo- 
mentum. 

Figure |lj demonstrates the excellent particle iden- 
tification for anti- protons, kaons, and pions within the 
50 GeV hadron beam using the Cerenkov counter under 
nitrogen, at a pressure of 160 torr. While the small sig- 
nal in the pedestal region is dominated by anti-protons, it 
could also be contaminated by other particles due to inef- 
ficiencies in C'l and C2 ■ An inefficiency study, performed 
by counting the number of pedestal events for the clean 
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Electron 


Hadron fraction(%) 


(GcV) 


fraction(%) 




5 


92 


8 


7.5 


72 


28 


15 


66 


34 (7r-:95.6, p+K--A.l) 


20 


< 1 


> 99 (7r-:95.5, p+K-:4.5) 


30 


< 0.25 


> 99.75 (7r-:94.9,p+K-;5.1) 


50 





100 (7r-:93.9, K-:3.1, p:3.0) 


75 





100 (7r-:91.7, K-:5.1, p:3.2) 


120 





100 (7r-:91, K-:6.2, p:2.8) 



Table 2: Summary of particle composition (e /tt /K /p) in the 
hadron calibration beam for various energies. 



muon sample with Cerenkov pressure above the muon 
threshold, shows that the Ci and C2 inefficiencies are 
less than 0.24% and 0.008%, respectively. 

Table |^ summarizes the beam particle composition 
for various hadron tunes. It is well demonstrated that 
the contamination from electrons in the beam for hadron 
tunes of momentum greater than 30 GeV is less than 
0.25%, minimizing the systematic error in hadron re- 
sponse measurements. 

8.4 The Calibration Spectrometer 

The spectrometer is designed to measure the absolute 
momentum of the calibration beam particles to better 
than 0.3% in an event-by-event basis. This is accom- 
plished by two means. First precisely calibrated dipole 
magnets are used, with J Bd£ known to better than 0.1% 
in the region traversed by the beam. Also, the bend an- 
gle is measured to better than 0.1% using drift chambers 
positioned over the 500 m beamline, which provides a 
long lever arm. This long length of the spectrometer 
chamber spacing allows us to tolerate a relative chamber 
alignment uncertainty of ^ 100 ^m. 

8.4-1 Upstream Tracking 

The position and the angle of the calibration beam tracks 
are determined at the upstream end of the spectrometer 
magnet by four 12 cm x 12 cm Single Wire Drift Cham- 
bers (SWDC)Ej. Each chamber consists of a pair of sense 
wires offset by ±2.03 cm from the beam center in each 
view. The operating gas, an equal mixture of Ar and 
C2Hg, and the field-shaping wires provide a saturated 
49 /im/ns drift speed over most of the gas volume. This 
can be seen in Figure ^ where the drift times measured 
on the two sense wires are plotted against each other. 
The dark band with slope — 1 is produced by tracks pass- 
ing between the two wires. The two bands with slope -1-1 
are produced by the tracks passing on the same side of 



-H 
IS 




Drift Time on x' Wire (nsec) 

Figure 15: T\ vs T2 for x— view of an SWDC plane. The cluster of 
points along the line with slope —1 to the left corresponds to the 
region between the x and x' wires; these events are used for track 
fitting. Nonlinear time-to-distance effects are visible for very long 
drift times and for events very close to the sense wires. 



the two wires. Small non-linear drift effects can be seen 
at very long and very short drift times. These effects are 
eliminated by using only tracks passing between the two 
wires. 

The chambers are grouped in two stations of two 
SWDC's each. One station is located immediately up- 
stream of the most upstream spectrometer magnet and 
the other is 83.3 m upstream of that station. 

Chamber position resolution can be determined by 
measuring the width of the distribution of differences in 
position measurements for a track passing between the 
two sense wires, and then dividing by ^/2. This is shown 
in Figure |l^, which demonstrates that the chambers have 
a resolution of 300 /im. Chamber alignment is achieved 
through an initial optical survey to an accuracy of 50 //m. 
The integrity of the alignment is constantly checked us- 
ing the straight test beam tracks, with the spectrometer 
magnets removed from their normal positions. Residual 
misalignments are estimated to be < 100 /xm and make 
a negligible contribution to slope and intercept measure- 
ments. 

Figure |l^ shows the distribution of slope, intercept, 
and for a typical upstream sample of tracks in each 
view. Tracks enter the spectrometer at projected angles 
of 9x = 74 mrad, 6y = —0.3 mrad with angular spreads 
of cTx = 0.2 mrad, ay = 0.08 mrad. The width of the 
beam is approximately 2.5 cm in x and y and is set by 
the trigger paddles. The distributions are consistent 
with their expected shape, giving us confidence that the 
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Figure 16: Difference in coordinates, Xi — X2 for an SWDC 
chamber plane from a sample of calibration beam tracks. The 
solid line represents a Gaussian fit on the distribution. The a of 
this distribution is 420 fim, implying a spatial resolution of 300 /im. 



upstream tracking is well understood. 
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Figure 17: A and B: Track slopes in x and y at the upstream 
end of caUbration spectrometer magnets as determined by cham- 
ber tracking. C and D: Beam profile in x and y at the upstream 
end of calibration spectrometer magnets as determined by chamber 
tracking. E and F: distribution of the upstream chamber fits in 
X and y views. Both distributions follow the expected normal 
distributions. 



8.4-2 Downstream Tracking 

The downstream section of the spectrometer consists of 
two 3 m X 3 m three-wire drift chambers separated by 
45.5 m. The first chamber downstream of the spec- 
trometer magnet string is positioned 23.7 m from the 
downstream end of the last momentum analyzing mag- 
net. Figure |l^ shows the distributions of the spec- 
trometer tracking fit, using the chambers upstream and 
downstream of the spectrometer dipole magnets. The 
fact that these distributions follow the expected dis- 
tribution gives us confidence in the absolute momentum 
determination. 

8.4-3 Spectrometer Magnet Calibration 

The Fermilab Majgnet Test Facility (MTF) calibrated the 
five EPB dipolesM (four plus an unused spare) used in 
the spectrometer. Precise / Bdi data are taken at the 
centerline of the magnet and are tied to magnet cur- 
rent and Hall probe voltage readout recordings. Shape 
studies are performed for J Bdl vs horizontal position 
at fixed vestical position and magnet shunt current mea- 
surementsta. Data are summarized by polynomial fits to 
the / Bdl measurements as functions of both Hall probe 
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Figure 18: distribution for X (A) and Y (B) view of calibration 
beam spectrometer track fit using chambers upstream and down- 
stream of the spectrometer magnet string. Ten points are used in 
the three parameter fit. 



12 



Yoke 



yoke 
iperature 



front view side view 

Figure 19: A schematic diagram of the Hall probe and the holder 
location inside the magnet. 



Magnet 


Probe 


Ai (m) 


Ao (T-m) 


A2 (T-'-m) 


11243 


95421 


-3.03617 


-0.00517 


0.00574 


11243 


95420 


-3.04009 


-0.00438 


0.00283 


11243 


95422 


-3.03576 


-0.00448 


0.00130 


11243 


95423 


-3.03392 


-0.00276 


0.00056 


11459 


95421 


-3.03770 


-0.00491 


0.00046 


11632 


95421 


-3.03575 


-0.00532 


0.00078 


11694 


95420 


-3.03983 


-0.00351 


0.00066 


20015 


95421 


-3.04207 


-0.00473 


0.00051 



Table 3: Coefficients of fits to J Bdl vs Hall probe readout for 
different calibration spectrometer EPB dipole/Hall probe combi- 
nations. The fits are of the form J BM = Ao + A^H + A2H^, 
where H is the Hall probe readout in Tesla and J Bdi is in units 
of Tesla-meters. 

output and magnet current. 

While shunt devices can be internaUy caHbrated to 
better than 1 part in 10^, the current reading in the 
two different configurations (i.e. different power sup- 
plies, buses, cables, and shunts) may differ by substan- 
tially more and cannot be used to obtain the absolute 
/ Bdi. Therefore, the absolute / Bdi is determined in 
the data taking configuration based on the Hall probe 
versus / Bdi calibration data. Further details of the 
magnet calibration are described elsewhere t^. 

Each of the four Hall probes is attached to a probe 
holder before the holders are mounted in the magnets. 
The probe holders are located approximately in the cen- 
ter magnet aperture on the lower pole face near the mag- 
net opening. This is shown schematically in Figure |l9|. 
Holders are angled to keep the probe cables from inter- 
fering with the beam and vice versa. Hall probes are 
read out and the values are recorded by the NuTeV data 
acquisition system once per spill, and J Bdi data are cal- 
culated for all events within a given spill using quadratic 
fits to the MTF Hall probe calibration data. 

Table ^ summarizes the J Bdi vs Hall probe fits 
to the data. The J Bdi-HaW probe relationship is very 
nearly linear, with the offset and quadratic corrections 
{Aq and A2 in Table ||) only a few parts-per-mil of the 
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Figure 20: Comparison of linear coefficients in J Bdi vs Hall probe 
readout for different EPB dipole/ Hall probe combinations for the 
NuTeV calibration beam spectrometer magnets. 



linear calibration constants Ai . Figure ^ compares the 
coefficient Ai for different probe-magnet combinations; 
variations from dipole-to-dipole and probe-to-probe are 
at the few tenths of a percent level. Linear fits do not 
parameterize the data to the required accuracy (< 0.1% 
deviation), but quadratic fits well describe both polari- 
ties of current ramping. Figure |2l| shows the fit results 
superimposed on the data for one of the magnets; plot- 
ted are the J Bdi points divided by the probe readout in 
order to accentuate non-linear effects. 

As a check, J Bdi values have also been calculated 
from high-order polynomial fits to J Bdi vs magnet shunt 
current data taken at MTF. Figure ^compares the Hall 
probe determination to the shunt current determination 
of three spectrometer magnets for a typical run. The two 
determinations agree within the expected precision of the 
shunt current measurement. 

8.5 The Calibration Trigger 

The calibration beam trigger consists of two small scintil- 
lator paddles shaped to shadow the "good field" regions 
of the spectrometer magnets. The two paddles are po- 
sitioned immediately upstream and downstream of the 
momentum analyzing magnet string. Figure ^ shows 
a schematic diagram of a calibration beam trigger scin- 
tillator paddle. The "good field" region, mapped out 
with J Bdi measurements described above, consists of 
the region across the face of the magnet over which the 
/ Bdi varies by less than 0.1% from its value at the cen- 
ter of the magnet. This unbiased trigger, with no energy 
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Figure 21: Results of quadratic fit to J Bdi vs Hall probe readout 
for a typical EPB dipole in the NuTeV calibration spectrometer. 
J BM/H vs H, where H is the Hall probe reading in Tesla and 

J BM in Tesla-meters is plotted to show the non-linear region at 
low fields. 
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Figure 22; Ratio of ^ BM calculated from magnet shunt current 
to the value calculated from the Hall probe for a typical NuTeV 
calibration run. 



9 Corrections for Systematic Effects 

In order to achieve sub-per cent precision in the abso- 
lute energy scale calibration, it is necessary for NuTeV 
to take into account a number of systematic effects. Most 
of these are due to the fact that the gains of the calorime- 
ter are determined using muons from the neutrino beam 
which averages the detector response over a long period 
of time (typically a week or more), while the calibration 
beam runs take place over much shorter time periods; for 
example, some hadron energy tests run for as little as an 
hour. The time dependent effects that need to be taken 
into account and that are discussed in this section are 
high voltage and temperature gain dependencies. 

There arc other differences between the energy de- 
position from a neutrino interaction and that from an 
incoming beam of particles. One of the differences is the 
calibration beam composition. The particle type depen- 
dence of the energy deposition is studied, and correction 
for the anti-proton contamination is applied for final en- 
ergy scale calibration. 

Another difference between neutrino and calibration 
beam energy deposition arises from the fact that a neu- 
trino may interact at different distances from a scintilla- 
tion counter, whereas the calibration beam always enters 
the calorimeter at the front and, in particular, electrons 
always interact in the first counter, which is preceded by 
two inches of steel. For this reason special care is taken 
to ensure that the calibration of the first few counters in 
the calorimeter is consistent with that of the latter ones. 
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9.1 Environmental and Voltage Monitoring System 

In order to obtain corrections for systematic effects, tem- 
perature, pressure, humidity, high voltages, and low volt- 
ages are monitored locally by a microprocessor controlled 
system. This system communicated periodically with the 
data acquisition computer to record the monitored data 
on the same tape as the neutrino and calibration data. 
The period for recording this data is one beam cycle, 
about one minute. 

The microprocessor signals were digitized by a 12-bit 
ADC and then read into a Basic Stamp BS2 microproces- 
sor. The signal is averaged over many readings to avoid 
noise. The results of this averaging were transmitted via 
a standard serial communication network to a personal 
computer in the control room. The personal computer 
monitored these values, issued warnings when values be- 
gan to deviate from the standard, issued alarms when the 
values were out of limits, and transmitted the raw data 
to the data acquisition system to be written to tape. 

The temperature of the calorimeter is monitored 
with a digital temperature integrated circuit to a res- 
olution of ±0.6°C. The temperature probes are placed 
in four locations on every fourteenth calorimeter counter 
unit. The absolute atmospheric pressure is monitored 
in four locations with a resolution of ±1.5%. The other 
environmental variables are monitored to a resolution of 
±2% in several places throughout the experimental hall. 

The calorimeter PMT high voltages are monitored 
using the LeCroy 1440 main frame readout system. The 
readout is recorded once per cycle to the data tape to- 
gether with neutrino and calibration data. 

9.2 Beam Component Correction 

The absolute hadron energy scale of the calorimeter is 
determined by measuring its response to single pion in- 
teractions. Any difference in the response of the detec- 
tor between pions and the kaons or anti-protons which 
contaminate the calibration beam must be accounted for 
when the absolute energy scale is set. These differences 
in the calorimeter response are investigated using clean 
samples of each particle type using the Cerenkov counter 
information. 

Based on the studies using the hadron beam at var- 
ious energies, we find that the calorimeter response to 
kaons agrees with that of pions. However, showers from 
anti-protons show higher responses (by ~1 GeV) than 
the showers from pions, due to the pp annihilation at the 
end of shower development process. This efifect has been 
discussed in a previous calorimeter review Ej. 

At high energies (> 50 GeV), the anti-proton effect is 
found to be negligible (< 0.03%), while this effect is very 
important at low energies. Table summarizes the sizes 



p 

(GeV) 


shower response 
(normalized to pion) 


fraction 
(%) 


correction 

to Ehad (%) 


5 


p 1.2 


3±1 


-0.6±0.2 


10 


p 1.1 


3±1 


-0.3±0.1 


15 


p+K--. 1.054 ±0.017 


4.1 


-0.22 


20 


p+K--. 1.033 ±0.010 


4.5 


-0.15 


30 


p+K": 1.027 ± 0.006 


5.1 


-0.14 


50 


p: 1.011 ±0.006 
K": 0.995 ±0.006 


3.0 
3.1 


-0.017 


75 


p: 1.008 ±0.004 
K": 0.997 ±0.004 


3.2 
5.1 


-0.010 


120 


p: 1.005 ±0.004 
0.998 ±0.003 


2.8 
6.2 


-0.002 



Table 4: Hadronic shower responses from kaon and anti-proton nor- 
malized to that of pions, and the correction factors to the hadronic 
shower energy, especially due to the anti-proton effect. 



of the correction factors, the contamination of kaons, 
anti-protons, and the shower responses (normalized to 
the pion shower response) of the calorimeter. 

9.S Muon Radiative Equilibrium (RE) Correction 

When a muon traverses material, it loses energy via 
electromagnetic processes: knock out electrons ((5-ray) 
from atoms, bremsstrahlung, e'^e~ pair production, etc. 
While most of the knock-on electrons are low energy elec- 
trons that do not penetrate deep into the material, high 
energy electrons from muon energy loss processes can 
leave energy in several counters. Thus, the energy de- 
posited in the most upstream few counters in the NuTeV 
calorimeter is relatively lower than other downstream 
counters since they have less material in before them. 
This effect is called the radiative equilibrium (RE) ef- 
fect. Since the NuTeV calibration beam enters the de- 
tector striking the most upstream counters, and the gain 
corrections for the counters are determined relative to 
muon energy deposit in a given counter, it is necessary 
for analyses to apply corrections to the gain factors to 
account for the RE effect. This effect causes an artificial 
over-estimate of the energy deposited in a few upstream 
counters relative to the downstream ones. 

We determine the size of this correction using a high 
statistics GEANT Monte Carlo study. Since this effect 
reduces the muon energy deposit gain normalization fac- 
tors for the most upstream and the second most upstream 
counters by 1% and 0.4%, respectively, the normalized 
hadron energy deposit in these two counters need to be 
reduced by the same factors. The resulting overall size 
of this correction to the hadronic response is less than 
0.1%. 
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9-4 Temperature Correction 

Many characteristics of the NuTeV calorimeter - PMT 
high voltage, PMT quantum efficiencies, scintillator light 
yield, electronics noise, etc - change with temperature. 
These changes contribute to the temperature dependence 
of the overall gain. Separating systematic effects from 
these different sources is difficult and unnecessary. The 
net effect that the temperature has on the overall gains of 
the PMT's can be measured by muon response maps and 
the average temperatures for each of these muon maps. 
That temperature dependence is shown in Figure || in 
Section ^. 

The muon map of the counter for a given running 
period is the time-integrated, beam-weighted response of 
the counter for a particular period of time. Since the 
neutrino data itself is, by definition, beam-weighted, the 
"average temperature" for a muon map is also the aver- 
age temperature for the neutrino data. 

However, this is not true for calibration beam data. 
A particular calibration beam study might only take a 
few hours, while the average muon map is calculated 
over a few weeks. If, during the few hours of the cali- 
bration beam study, the temperature is significantly dif- 
ferent from the average temperature for the muon map 
used, the gain of the calorimeter during that brief time 
interval would be different from the muon map average 
gain. To correct for these gain differences, a tempera- 
ture correction is applied to the calibration beam data 
on an event-by-event basis, such that the effective muon 
map used would be the appropriate muon map for that 
particular temperature. 

To calculate the corrections, the average gain for each 
counter (G = 9^^+3B+gc+gD ^ ij^garly correlated to the 
measured temperature from the sensors located on the 
calorimeter (G(T) = A x T + B). Although there are 
sensors placed along the length of the calorimeter and 
the temperatures are measured throughout the experi- 
mental hall, the temperatures measured by the sensors 
near the least insulated part of the calorimeter are used to 
determine the temperature correction. The temperature 
correction to each counter is simply the first term in the 
Taylor expansion of ratio of G{Tmuon map)/G(Tcurrent), 
or (1- § X {Tmuon map - Tcurrent)) ■ Temperature correc- 
tions tend to be as large as several tenths of a percent. 



9.5 High Voltage Correction 

Counters in the NuTeV calorimeter have four PMT's, one 
in each corner, as described in Section |^. The overall gain 
of a given counter depends strongly on the combination of 
the individual PMT gain. One of the systematic factors 
that directly affects the gain is PMT high voltage (HV). 
Thus, for a high precision calibration, it is important 
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Figure 24: Typical HV readout values for four randomly selected 
PMT's. As one can observe, the HV readout does not vary more 
than 2-5 V. 



to correct for overall gain fluctuations due to any HV 
variation. 

The PMT gain variation as a function of HV is mea- 
sured prior to running for all PMT's used in the calorime- 
ter, and is parameterized as: 



gpMT — aV"', 



(8) 



where gpMT is the gain of the given PMT, V is the HV 
in units of volts, and a and a are the fit parameters. The 
exponents a are determined for each PMT, and a typical 
value of a is ~ 6.8. 

The NuTeV experiment implemented a slow mon- 
itoring system th at m onitored PMT HV values as de- 
scribed in Section 9.1. Six LeCroy 1440 HV mainframes 
supply high voltage to the calorimeter PMT's. The slow 
monitoring system read out 1 HV channel per second and 
completely cycled through all HV channels in 5 minutes. 
The entire record in the database is written to the neu- 
trino data tape once every beam cycle as the last record 
in the given cycle. The readout resolution of the NuTeV 
HV slow monitoring system is ~ 1 V; the PMT's are 
typically set at between -1400 V and -1500 V. 

A study based on a total of 280,000 measurements of 
each individual HV read-back, taken over the entire run 
period, reveals that the typical variation of each PMT 
HV readout is within 2 V and the RMS of the distribution 
is typically less than 0.5 V. Figure p4 shows typical HV 
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readout values of four randomly selected PMT's for all 
the readings throughout the entire run. 

Despite the fact that we expect very small corrections 
due to HV variations based on the HV readout measure- 
ments, we correct for HV for the calibration beam due 
to the fact that calibration runs are typically localized in 
time while the gain correction factors are averaged over 
a longer time period. 

HV corrections for each counter are done using the 
measured parameters in Eq. ||, relative to the counter 
gain correction factor at the center of the counter, av- 
eraged over the muon response map period. The rela- 
tive gain correction factor is computed using Eq. ^ on 
an event-by-event basis and normalizing the gain at the 
given HV readout to that at the average HV readout of 
the PMT in the given run. The relative correction factor 
for counter z, fl^y, is defined as: 

jhv - — ' (y) 

where < Vij > is the average HV readout value of PMT 
J of the counter i for the given run, Vij is the HV readout 
for the given event, aj is the exponent in Eq. ^ for PMT 
j, and PHij is the individual pulse height from the PMT 
J. As we expect, the typical overall size of this relative 
correction factor is on the order of 0.1% or less. 

10 Measurement of Muon Energy Loss in the 
Calorimeter and Comparison with GEANT 

The toroid spectrometer is located downstream of the 
690-ton NuTeV calorimeter. For an accurate measure- 
ment of each muon's momentum, the energy lost by the 
muon in the calorimeter (A£') has to be included. A 
precise measurement of AE is also necessary for the cal- 
ibration of the toroid using test beam muons. Knowl- 
edge of muon energy depositions is also needed for the 
hadronic energy measurement, since muons, originating 
in charged current interactions, contribute to hadronic 
shower pulse heights. 

A minimum-ionizing particle passing through the de- 
tector leaves a characteristic energy deposit in each of 
the scintillation counters. The energy loss of a muon 
traversing the calorimeter changes with energy. For high 
energy muons the contribution to the muon energy loss 
from bremsstrahlung, electron-positron pair production, 
and nuclear interactions increases. These processes may 
yield larger electromagnetic showers than would be true 
for a strictly "minimum-ionizing" particle. A coarse sam- 
pling calorimeter, such as NuTeV's, is strongly non-linear 
in energy between a fraction of an MeV to a few GeV. 
Hence the conversion of the light yield induced by a muon 
passing through the counters to measured muon energy 



Process 


Value 


Rayleigh Scattering 
7-Induced Fission 
(5-ray Generation above 
Restricted Landau below 
Direct Pair Production 
Bremsstrahlung Tracking 

Other Particles 


ON (IRAYL=1) 
ON (PFIS=1) 
DCUTE = 100 keV 
DCUTE = 100 keV 
ON (PPCUTM = 2.04 MeV) 
(BCUTE = 100 keV) 
(BCUTM = 100 keV) 
CUTGAM^lOO keV 
CUTELE = 100 keV 
CUTNEU = 100 keV 
CUTHAD = 100 keV 
CUTMUO = 1 MeV 



Table 5: Parameters with changed values from their default values 
in GEANT V3.215. 



loss in the calorimeter requires a differentiation between 
lower and higher energy processes. We use a GEANT- 
based simulation of the detector (McNuTeV) to deter- 
mine the best pulse height to GeV conversion method 
or "reconstruction algorithm" for the total energy lost 
by the muon, AE. Tests of both the GEANT simula- 
tion and the reconstruction algorithm are described in 
the following sections. 

10.1 Counter Pulse Height Simulation using GEANT 

The NuTeV calorimeter simulation segments the 
calorimeter into six identical carts, each of which con- 
sists of seven unit calorimeter layers described in Sec- 
tion ^ The steel, water, drift chamber gas. Incite, mylar, 
polythene, air, copper, and GIO are specified as separate 
GEANT volumes building the layers with sizes and con- 
figurations closely matching the physical detector. We 
find that very detailed modeling of the detector is neces- 
sary to achieve good agreement between calibration beam 
data and the GEANT simulation of muon responses. 

We use version 3.215 of GEANT and set the physics 
control variables to their default values, with the excep- 
tions listed in Table ||. The energy deposited in scintilla- 
tion counters follows an avalanche in our 10-stage model 
of PMT's and is smeared statistically at each step. The 
number of photoelectrons used in the smearing is tuned 
to match the widths of muon dE/dx deposition in the 
data. Pedestals, gains, and the digitization of LOW and 
HIGH channels of electronics are also simulated. Ob- 
served pulse heights, both for data and for the simula- 
tion, are expressed in units of MIP's, where 1 MIP is 
defined as a truncated mean of the energy loss of 77 GeV 
muon (see Section ^ . The resulting GEANT events are 
passed through the same analysis chain as the actual data 
events. 
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Figure 25: Counter response to muons in calibration beam (his- 
togram) and Monte Carlo (crosses) for muon of energies of 15 (A), 
50 (B), 100 (C), and 166 GeV (D). DVTEN (HIGH/LOW channel 
depending on saturation) are muon map corrected and measured 
in MIP's. All 84 counters contribute to the distributions. Solid 
lines represent 5-parameter asymmetric Gaussian fits (Eq. |lo|) to 
the distributions. 



10.2 Data/GEANT Comparisons for Muons 

Muon calibration beam data are taken throughout the 
1996-97 NuTeV run totahng approximately 250 10000- 
event data sets. Most of that data are 50 GeV muon sets 
used for measuring the magnetic field of the toroid. An- 
other subset, also used in this study, consists of runs with 
muon energies spanning from 12.5 to 190 GeV. GEANT 
samples are generated with the energies, momenta, and 
positions at the entrance to the calorimeter matching the 
calibration beam data samples. In the comparisons, cuts 
are applied to the calibration beam muons to assure that 
the momentum measured in the test beam spectrome- 
ter and the x- and y-vertex positions are reconstructed 
within ±3 standard deviations around the mean value. 

Figure ^ illustrates the detector response to cali- 
bration beam data (histogram) for muons of energies of 
15, 50, 100, and 166 GeV. The GEANT simulation is 
marked by crosses. DVTEN shown in the plots is the 
pulse height of a counter measured in units of MIP's, 
after application of the position dependent gain correc- 

It uses the HIGH chan- 



tion discussed in Section 6.1 



nel of electronics until its ADC saturation (1900 ADC 
counts), and the LOW channel readout above that. The 
DVTEN distributions in muon energy bins are fitted with 



Figure 26: Fit parameters (A, B: mean and RMS of the Gaussian 
part, C and D: asymmetric tail, Eq. hoh to DVTEN distributions 
for Monte Carlo (open circles) and caliBration beam muons plotted 
versus incident muon energy (E^). 



a five-parameter asymmetric Gaussian fit, !F, in which 
the width of the Gaussian runs on one side, varying with 
the X-axis: 



P3max(l, (1 



{Pi + xPr,){x - P2))) 



(10) 
(11) 



Figure |2^ gives the values of the four parameters of 
these fits (peak of the Gaussian portion of the distribu- 
tion P2, its width P3, and, P4 and P5, two parameters 
describing the asymmetric tail) as a function of muon 
energy. The Monte Carlo is represented by open circles, 
while the muon data sets by solid circles and stars. The 
calibration beam data show that the most probable value 
of muon energy loss in a counter is independent of muon 
energy in the 20-200 GeV range (see Figure ^6|) . This is 
to be contrasted with almost linear increase with muon 
energy of the mean energy loss in the counter (also Fig- 
ure p9|) . Figure ^ shows the pulse heights summed over 
84 counters traversed by muons (DVTEN SUM distri- 
butions). This summed plot would magnify small dis- 
crepancies (e.g., in the tails) in comparison of calibration 
data and the Monte Carlo, but GEANT still describes the 
data well. A summary plot containing fit parameters to 
DVTEN SUM histograms for all available muon energy 
points is shown in Figure Eq. The mean and RMS val- 
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Figure 27: Counter response to muons summed over 84 counters in 
calibration beam (histogram) and Monte Carlo (crosses) for muon 
of energies of 15 (A), 50 (B), 100 (C), and 166 GeV (D). 
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Figure 28: Fit parameters (A, B: mean and RMS of the Gaussian 
part, C and D: asymmetric tail, Eq. |lo| to DVTEN SUM distribu- 
tions for Monte Carlo (open circles) and calibration beam muons 
plotted versus incident muon energy (E^). DVTEN SUM is the 
sum of the DVTEN over 84 counters. 



ues (not from fits, but from the histogram statistics) for 
the DVTEN and DVTEN SUM distributions are plotted 
in Figure From these caHbration beam data-GEANT 
comparisons, we conclude that both the low and high en- 
ergy components of the muon energy loss in the calorime- 
ter are well modeled in our Monte Carlo over the full 
scale of such depositions. Modeling of GEANT muon 
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Figure 29: Mean values of DVTEN (A) and DVTEN SUM over 
84 counters (C) distributions for Monte Carlo (open circles) and 
calibration beam muons (solid circles) as well as RMS values of 
these distributions. 

energy loss in the steel and remaining absorber materials 
can be checked using so-called range-out muons. Those 
are low energy calibration beam muons that stop in the 
calorimeter. Our measurements give the mean length of 
75.7±0.3 counters for a 12.5 GeV muon in the data ver- 
sus 74.7±0.09 in the simulation. The agreement between 
data and Monte Carlo in the mean length is only within 
about three standard deviations. One possible cause of 
this disagreement is the lack of any correlated counter 
noise in the Monte Carlo. 

10.3 Study of Muon Energy Loss Reconstruction Algo- 
rithm Using GEANT 

The goal of this study is to find an optimal algorithm 
to determine the muon energy loss in the calorimeter, 
Ai?, from the observed muon pulse height. We recon- 
struct Ai? and compare it to the "true" Ai? known from 
GEANT on an event-by-event basis in the range of 15- 
190 GeV. We define TOUT as the counter pulse height 
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Figure 30: Method a: The width (A) and mean (B) from Gaussian 
fit (sohd circle) and histogram statistics (stars) of the difference 
between "true" AE and reconstructed AE distributions in muon 
energy loss bins. C and D: the same in muon momentum bins. 



at which we switch from applying ( "low energy" ) to 
Ce ("high energy") conversion from MIP's to GeV, un- 
der the assumption that sufficiently high pulse heights 
arise from electromagnetic processes sampled over sev- 
eral counters. Three different reconstruction schemes for 
AE are studied: 

Method a: A model of two conversion constants 

and Ce and TCUT of 3 MIP's (this is a scheme 
used in our predecessor experiment CCFR, where 
Ce is determined from electron calibration beam 
data) . 

Method (3: A "one-function model", where one func- 
tion Cfj_{E^) is used to account for ionization and 
the increase of the radiative component of dE/dx 
with energy. 

Method 7: A model of conversion function C^{E^) 
varying with muon energy i?^, applied below 
TCUT of 5 MIP's, and a constant conversion Ce 
above that TCUT. 

As an illustration, we show the widths and the means of 
the difference between the "true" and reconstructed Ai? 
distributions in either muon energy bins or in the bins of 
the muon true energy loss for Method a in Figure ^ and 
Method 7 in Figure Both the fitted mean and sigmas 
of Gaussian fits to these distributions (solid circles) and 
the average values and RMS (stars) of these distributions 
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Figure 31: Method 7. The fitted widths (A) and means (B) of the 
difference between true AE and reconstructed AE distributions in 
muon energy loss bins. (C, D): the same in muon momentum bins. 



are plotted. Notice that Method a underestimates AE 
for high energy electromagnetic depositions if the lower 
end of the energy loss spectrum is set to match GEANT's 
"true" AE (Figure |^(B)). Similarly for Method (i (not 
shown) - no C^[E^) can be found that describes the 
conversion from MIP's to GeV for both the most prob- 
able and the mean dE/dx at the low and high ends of 
muon energy spectrum at the same time. In Method 
7, where the variation in low-energy radiative deposi- 
tions with muon energy is accounted for by variation of 
the conversion function C^(i?^), we find the best match 
of the "true" AE for all muon energies and AE values 
(Figure |3ll). The function C^(i?^), based on the best cal- 
ibration beam-to-GEANT match, is shown in Figure 

Figure ^ illustrates the total muon energy loss over 
the length of NuTeV calorimeter in terms of mean and 
the most probable value, where the latter is defined as 
the result of a fit of asymmetric Gaussian function 
(parameter P2) to the AE distribution. Figure ^ gives 
a ratio of the most probable AE for 50 GeV calibra- 
tion beam muons, traversing the NuTeV calorimeter at 
different angles and transverse positions, to a nominal 
15.2 GeV GEANT prediction for their energy loss. The 
ratio is plotted as a function of muon azimuthal angle, (f), 
at the most upstream surface of the detector. As can be 
seen in this figure, we reconstruct muon AE to within 
±0.7%, independent of 4> (or position in the counter). 
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Figure 32: Energy dependence of calibration beam muon deposi- 
tions of pulse heights below 5 MIP's, over 84 counters, shown in 
units of GeV/MIP and the parameterization C^(i?fi). 



This is an important verification of the muon map cor- 
rection and counter gain stability over time. The counter 
pulse heights for these calibration beam muon samples 
are corrected by gains that vary by as much as an order 
of magnitude, depending on their transverse vertex and 
pathway through the calorimeter. 

11 Shower Energy Definition 

The definition of the shower energy used in the hadron 
energy calibration is: 



E. 



shower 



place place-{-19 
i—1 i— place 



, (12) 



where i is the counter number; PHi is the pulse height 
normalized for the muon gain as described in Section |6| in 
units of MIP in counter i; hi is the hadron/muon gain ra- 
tio (described in Section Q); and Ctt is the hadron calibra- 
tion constant. The variable pZace is the counter where the 
hadron shower started to develop, and is determined by 
an algorithm designed to locate where a neutrino interac- 
tion begins, place is defined as the upstream of two con- 
secutive counters which have more than a certain number 
of MIP's, where that number depends on the total energy 
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Figure 33: Most probable (A) and mean value (B) of total muon 
energy loss in NuTeV calorimeter. Comparison of GEANT predic- 
tion (open circle) and calibration beam measurement (solid circles) 
versus muon energy. 
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Figure 34: Ratio of the most probable values of AE for 50 GeV 
calibration beam muons to GEANT prediction (15.2 GeV) as a 
function of muon azimuthal angle. 



of the hadron shower and is at least four. Upstream of 
place the hadron is treated as a minimum ionizing par- 
ticle. In contrast, electron showers always start at the 
first counter. For the electron energy measurements de- 
scribed later only the seven most upstream counters are 
used, and the hadron/muon gain ratio is applied to each 
counter's pulse height. 

In order to determine the hadron energy calibration 
constant which, in our definition, is the GeV-to-MIP con- 
version factor (C^), it is important to define the energy 
variable to contain the entire shower. On the contrary, 
for a precision measurement, one does not want to in- 
clude too many counters in the sum because adding more 
counters than necessary would introduce noise into the 
system due to pedestal fluctuations in the counters. Fi- 
nally, the algorithm should be as close as possible to that 
used in analyzing the neutrino data, which also sums over 
a certain number of counters following p/ace. We perform 
a study to optimize the number of counters over which to 
sum the pulse heights. Using hadron beams over the en- 
ergy range between 10 GeV and 190 GeV, we determine 
that summing over 20 counters beginning from the most 
upstream counter is optimal for calibration purposes. 
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Figure 35 shows the cumulative fractional energy as 
a function of the number of summed counters for vari- 
ous hadron beam energies. Since we summed the pulse 
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Figure 35: Cumulative fractional shower energy as a function of 
the added number of counters. Note that for all energies, adding 
20 counters provides full longitudinal shower containment. 



heights of at least 20 counters dependent on the mea- 
sured hadron energy, we introduce an energy dependent 
noise level. This noise level depends on beam energy 
because the shower penetration depth depends on beam 
energy. The number of counters without actual shower 
energy increases with decreasing beam energy and the 
noise level gets more prominent for low energy beams. 
Therefore, the low energy calibration has a larger con- 
tribution from this noise effect. However, since hadronic 
energy resolution is worse at low energies, this noise is 
less important. 

12 Hadron Energy Response and Resolution 

The simplest test of the muon calibration technique de- 
scribed above is the time dependence of a particular cal- 
ibration beam setting. Figure ^ shows the time depen- 
dence of both 50 GeV hadron and 100 GeV hadron runs 
that are taken periodically during the course of the ex- 
periment. The RMS of the ratio between reconstructed 
calorimeter energy and beam momentum is 0.4-0.5%, and 
is due to the statistical uncertainty in the muon maps 
themselves, as well as the electronics gain coefficients. 

The calorimeter response to a monochromatic beam 
of hadrons can be characterized by a function similar to 
a Poisson distribution. This is because the energy recon- 
structed by the calorimeter is proportional to the num- 
ber of shower particles produced by the incident hadron. 
The statistical fluctuation of the number of electromag- 
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Figure 36: Remaining time dependence of the ratio of reconstructed 
calorimeter energy divided by measured beam momentum. The left 
graph is for 50 GeV runs, and the right is for 100 GeV runs. The 
horizontal axis is time in units of blocks, where each run block cor- 
responds to a period of about 2 weeks. The responses are obtained 
after all the time-dependent corrections are applied to the data, 
but before the final energy scale is set. 



netic particles in the shower causes the response to look 
Poisson-like at low energies, and to become Gaussian at 
high energies. This can be seen in figure 37, where the 
5 GeV data are much less symmetric than the 190 GeV 
data around the peak of the Energy/momentum distri- 
bution. Figure ^ shows that the Poisson-like function, 
defined below, describes the data over several decades. 
Fluctuations in where the primary hadron interaction oc- 
curs can also contribute to the asymmetric shape of this 
distribution, but again these fluctuations have a negligi- 
ble effect at high energies. 

The Poisson distribution is normally written as 

PiN,^i)^^^——, (13) 



iV! ' 

where P{N, ii) is the probability of seeing N shower par- 
ticles if fj, are expected. The RMS of this distribution is 
^/JI and the mean is /i. As iV gets large this approaches a 
simple Gaussian distribution. To remove any effects from 
variations in time of the calibration beam momentum, 
the calorimeter energy is divided by the reconstructed 
particle momentum on an event-by-event basis. This im- 
plies that the mean of the distribution is decoupled from 
the width, but the fractional width (width divided by the 
mean) remains 1/^/1. Generalizing Eq. |l^ to decouple 
the mean from the width and expanding about the peak, 
we can parameterize the Poisson distribution as follows 
(note that keeping only the first term in F{x) gives a 
Gaussian distribution) : 

-F{x) 



where 



and 



P{x) = Ae' 



(14) 



(15) 
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Figure 37: Poisson fits to the calorimeter energy divided by mo- 
mentum distributions for four different energies: 5 GeV, 10 GeV, 
75 GeV, and 190 GeV. 
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where E is the measured hadron energy, p is the recon- 
structed cahbration beam momentum, the peak E/p — 1 
of the distribution is C, and the width of the distribution 
is B. The fractional width of the distribution is 1/a/B. 
At beam momenta of 30 GeV and above, this equation 
is very close to a Gaussian distribution. Figure 37 shows 
fits to the above equation for four different energies. 

The energy dependencies of the mean E/p distribu- 
tion and the Poisson widths are shown in Figure 39. If 
the hadron calibration constant C^, as defined in Sec- 
tion is set to 0.212 after all corrections, the mean 
energy response divided by the reconstructed test beam 
momentum at 75 GeV is l.OOOiO.OOl. Note that the non- 
linearity of the calorimeter between 10 GeV and 190 GeV 
is only 3%. This comes from the fact that electrons and 
hadrons have a very similar response, and so the electro- 
magnetic component of the shower, which changes as a 
function of energy, will not change the reconstructed en- 
ergy. The Poisson widths can be fitted to the standard 
form a{E)/E = A(B ® ^, where A is a constant term 
coming from calibration uncertainties, B is the stochas- 
tic term from the sampling of the shower, and C is from 
noise due to pedestal fluctuations. The data show no ev- 
idence for a noise term and so C is removed from the fit. 
The stochastic term is proportional to the square of the 
thickness of the sampling layer. 

The energy points below 10 GeV have to be mea- 
sured using a different energy algorithm, since the one 
that is designed for neutrino interactions (requiring at 
least 4 MIP's in two consecutive counters) is biased for 
hadron showers below 10 GeV. For these lowest points, 
the energy is determined simply by summing the most 
upstream twenty counters, and is using the hadron gain 
coefficients for each counter. To remove electrons from 
the low energy samples, cuts are made based on Gerenkov 
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Figure 39: Hadron energy response versus reconstructed test beam 
momentum and comparison with fit to Groom's parameterization 
for non-Hnearity, and Poisson width distribution versus energy with 
fit to a{E)/E = AS>^. The open symbols are lower energy runs 
with slightly difi'erent cuts and are not used in the fits. 
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Figure 38: The distributions of calorimeter responses to 50 and 
100 GeV hadrons with very high statistics. The data follow the 
Poisson-like shape over several decades. 



24 



E 1.02 




e/h=1.0^ 
e/h=1.2 



in the hadron/muon gain ratio for the counters. 



5 6 7 8 9 10 20 
Momentum (GeV) 

Figure 40: Lowest energy hadron response versus reconstructed 
test beam momentum and comparison with fit to Groom's param- 
eterization for non-linearity. The open symbols axe lower energy 
runs with slightly different cuts. 



counter particle identification system in the beamUne. 
To remove muons in the hadron beam, a loose cut is 
made on the most upstream of three consecutive coun- 
ters that have less than 0.25 MIP's in them. The lat- 
ter cut removes events caused by muons in the hadron 
beam, but did not remove events with secondary muons 
created in the hadron shower. Finally, because the low- 
est energy points have low statistics the means of the 
energy /momentum distributions are plotted rather than 
the results of the Poisson fits. Figure pO| shows the nonlin- 
earity of the NuTeV calorimeter to low energy hadrons. 
For energies 5.9 GeV and above. Groom's parameteriza- 
tion (see section |l^) with e/h— 1.08 (solid line) agrees 
well with the data. The overall agreement with the pa- 
rameterization is not improved by changing e/h from 1.08 
in the hadron response curve, as is also shown in figure 
1- 

Finally, any additional position dependence not 
taken into account in the muon map procedure outlined 
earlier is studied using a 75 GeV hadron beam aimed 
at diff'erent locations on the front face of the calorimeter 
using the rotating dipole at the end of the momentum 
analyzing magnet chain in the spectrometer. For hadron 
showers that start more than 50 cm from the closest edge 
of the detector, the energy reconstruction is constant to 
better than a 0.5% in the calorimeter response, when 
normalized using the muon maps. By aiming the hadron 
beam as close to the edges as is safe, it is determined 
that hadron shower leakage does not begin to affect en- 
ergy reconstruction until the shower starts at 25 cm from 
the calorimeter edge. 

Table || lists systematic errors that contribute to 
the uncertainty in overall hadron energy scale of the 
calorimeter. It is clear from the list that the largest sin- 
gle systematic error is due to the statistical uncertainty 





Fractional 


Source 


Uncertainty 


Hall Probe Readout 


0.03% 


(from shunt comparisons, see Figure ^) 




Magnetic Field Homogeneity 


0.03% 


(from Position-dependence measurements) 




Beam Composition Corrections 


0.03% 


(100% of effect above 30 GeV) 




Transition Effect Uncertainties 


0.03% 


(10% of effect on hadrons) 




Temperature Corrections 
(10% of effect) 


0.02% 


High Voltage Correction 


< 0.01% 


Spectrometer Alignment 


0.1% 


Uncertainty in hadron/muon gain ratios 


0.4% 


Fit Normalization Error 


0.1% 


Statistical 


< 0.01% 


Total 


0.43% 



Table 6: Table of contributions to uncertainty in overall hadron 
energy scale. 



13 Hadron Response Comparison to Monte 
Carlo 

The task of reproducing the calorimeter attributes in a 
GEANT-based Monte Carlo is a challenging one. This is 
demonstrated in Figure Both the non-linearity and 
hadron energy resolution of the calorimeter depend crit- 
ically on the difference in the calorimeter's response to 
hadrons and to electrons (c/tt). In order for a Monte 
Carlo to simulate hadrons correctly, it must first cor- 
rectly simulate the calorimeter's electron response, and 
then have an accurate hadron shower model, a thorough 
description of the geometry of the calorimeter, and an ac- 
curate model for the way particles propagate in the par 



ticular media that comprise the calorimeter. Section 14 
describes the calorimeter electron energy response and 
the resolution in detail; the studies show that the electron 
energy resolution is well-modeled in the detailed GEANT 
simulation of the calorimeter. 

Once the electron to hadron response is measured 
in the data at a particular energy, one can minimize 
the dependence on hadron shower generators by only 
using them to predict the fraction of tt'^'s produced in 
a hadron shower as a function of energy, f^a{E). Fig- 
ure ^ shows three different hadron generators' predic- 
tions for the fraction of electromagnetic energy deposited 
in a hadronic shower as a function of energy. Two param- 
eterizations for tbia fraction are also shown; Wigmans' 
parameterization E2l is /^o = 0.11 ln(£') and Groom's is 
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Wigmans:f|^„=0.11ln(E) 
Groom: f -1.-(E/0.96)'"'^'* 
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Figure 41: Fraction of electromagnetic energy in a hadronic shower 
as a function of hadron energy for three different hadron shower 
generators: GHEISHA, GFLIWA, and GCAIrQR. Fits of the sim- 
ulation to both the Wigmans'E3 and Groom'stZI parameterizations 
are also shown. 
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Figure 42: Predicted non-linearity for three different hadron 
shower generators: GHEISHA, GFLUKA, and GCALOR, where all 
three are required to have the same value of E / E{true) at 50 GeV. 
The prediction for both the Wigmans and Groom parameteriza- 
tions are also shown. 



The reconstructed energy of a shower is defined as 

E = EtrueiHl - Uo{E)) + eUo{E)), (17) 

where h is the ratio of reconstructed to true energy for 
a "pure" hadron, and e is that same ratio for a "pure" 
electron. The non-Hnearity as a function of energy can be 
expressed as the ratio of R{E) / R{Eref) , where R{E) = 
E/Etrue- In other words, 

non - hnearity = i_f;^^(-°|^'^]+|[;°gL.) ■ (18) 

By requiring the three generators shown in Figure ^ 
to have the same R{E) / R{Eref) at 50 GeV, one can con- 
struct the expected non-linearity as a function of energy. 
This is shown in Figure Although the fraction of elec- 
tromagnetic energy at a given energy varies among the 
different generators, the predicted non-linearity is simi- 
lar. 

By fitting the hadron energy response shown in 
Figure |3^, assuming the non-linearity predicted by 
Groom's parameterization, one arrives at a value of of 
1.079 ± 0.011. Given that at 75 GeV the fraction of t^^'s 
is roughly 50%, f - 1 = 0.5(f - 1). Therefore, the elec- 
tron to hadron response ratio is roughly of 1.035 ± 0.01 
at 75 GeV. This is in agreement with what is seen in 
the calibration beam comparisons between hadrons and 
electrons, as will be described in Section O. A similar 



conclusion has been reported by the CDF coUaborationEj 
for their plug-upgrade hadron calorimeter, which has a 
much larger non-linearity. 

As a side note, if we use the different hadron shower 
generators to predict the ratio of electron to hadron re- 
sponse at 75 GeV, they all give different ratios. Since 
the difference between these responses is also a large fac- 
tor in determining the hadron energy resolution, they all 
predict correspondingly different hadron energy resolu- 
tions (the closer the electron to hadron response ratio is 
to unity, the better the hadron energy resolution). 

GFLUKA E3 predicts the lowest ratio of electron to 
hadron response, and also predicts a hadron energy res- 
olution of Qm/^E{GeV). GHEISHAB predicts a ratio 
of electron to hadron response of 1.15 at 75 GeV rather 
than 1.09, and predicts a hadron eofirgy resolution of 
l.l^/^E{GeV). Finally, GCALOR eI predicts that the 
ratio of electron to hadron response is less than unity, 
rather than greater than unity. 

All three generators are tested using identical 
GEANT energy cutoff settings and identical geometry 
input. Although GHEISHA is native to the GEANT pro- 
gram, GFLUKA and GCALOR are imperfect implemen- 
tations of the original FLUKA and CALOR programs 
and have been known to produce jSomewhat different re- 
sults than the original generatorsEil. 

Finally, although the hadron non-linearity is now 
parameterized and well-measured, the purpose of the 



26 



NuTeV calorimeter is to measure hadron showers gen- 
erated by neutrino interactions, not hadron showers gen- 
erated by a single charged hadron. To study any possi- 
ble difference the LUND Monte Carlo program is used 
to determine the first set of particles produced from the 
hadron shower of a neutrino interaction. Groom's param- 
eterization is then used to calculate the electromagnetic 
fraction of the charged hadrons which get produced. Al- 
though the charged hadrons have lower energy than the 
initial total hadron energy, and as such would have a 
lower electromagnetic fraction, there are also neutral pi- 
ous that are produced, which increase the electromag- 
netic fraction. The two effects cancel, keeping the elec- 
tromagnetic fraction as a function of total hadron en- 
ergy the same between neutrino-induced hadrons and sin- 
gle hadrons, to a few per cent. The resulting effect on 
the neutrino-generated hadron non-linearity is negligible 
compared to the statistical error on the e/h fraction it- 
self, and the overall energy scale change is consistent with 
zero to better than 0.1%. 

14 Electron Energy Response and Resolution 

The calorimeter response to electrons can be measured 
using the calibration beam when set to the electron 
mode, as described in Section ^.2[ . Although there is 
a large contamination of muons in the electron running, 
this is easily removed from the data sample by looking 
at the most downstream counter and by selecting events 
with more than one minimum ionizing particle in them. 
Since electrons penetrate no more than a few cm of steel, 
most of the energy is deposited in the first three coun- 
ters, so calibrating the detector response to electrons is 
extremely dependent on the muon maps of those three 
counters, and has larger systematic uncertainties due to 
statistical errors in the muon maps. 

Figure ^ shows the shape of the electron energy de- 
position for 30 and 170 GeV electrons, as well as the fit 
to the Poisson distribution. The detailed GEANT sim- 





36.71 / 39 
139.5 
327.7 
-0.8397E-02 

170 GeV 
e 



0.5 



1.5 2 
Energy/Momentum 



0.5 1 1.5 2 

Energy/Momentum 



Figure 43: Calorimeter energy divided by test beam reconstructed 
momentum distributions for 30 GeV and 170 GeV electrons, and 
the results from the Poisson fit to the distribution. 
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Figure 44: Normalized response of calorimeter to electrons as a 
function of energy for both data (solid circles) and GEANT-based 
Monte Carlo (open squares). 



ulation reproduces the shapes of the distributions well, 
but there are substantial differences between the mean 
values as a function of energy. 

Figure ^ shows the electron energy divided by mo- 
mentum for both data and Monte Carlo. There are sev- 
eral contributions to the large non-linearity. The most 
important contribution to the non-linearity, and the only 
one present in the GEANT simulation, is due to the 
fact that electrons arrive at the upstream edge of the 
calorimeter and begin showering immediately. If one gen- 
erates the electron showers in the simulation equally dis- 
tributed throughout the steel, as neutrino interactions 
(and hadrons to first order) would be, the GEANT re- 
sponse is linear to better than 0.5%. 

There are two other effects that are present in the 
data but not in the simulation, and these give rise to 
the additional non-linearity that is seen in the data. The 
electron response depends critically on the both the muon 
gain and the hadron/muon gain ratio for the first two 
counters. These ratios are only known to about one per 
cent, per counter, so this contributes an additional un- 
certainty which could affect both the scale and the non- 
linearity. Another contribution to the non-linearity is due 
to the fact that there are cuts on the reconstructed track 
momentum in the data that cannot be made in the Monte 
Carlo simulation. These cuts do not affect the hadron re- 
sponse because hadrons are extremely unlikely to shower 
before the calorimeter. There is however approximately 
one radiation length of material distributed throughout 
the last arm of the calibration beam spectrometer. This 
material is included in the GEANT simulation, but its ef- 
fect on the tracking efhciency is not. This is particularly 
important at high energies, where the upstream showers 
are most likely to cause ambiguities in the momentum 
determination. 

One correction to the gains of the first few counters 
in the detector that is extremely important for measuring 
the electron response of the calorim eter , is due to the RE 
effect which is discussed in Section 
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Figure 45: Electron energy resolution for both data and GEANT- 
based Monte Carlo as a function of energy, and the results to a fit 
of the form cj{E)/E = A® 



The GEANT-based Monte Carlo predicts the sam- 
pHng term in the resolution of electron energies to within 
2%, as can be seen in Figure |4^. As with the hadron res- 
olution, the electron resolution can be fit to the form 
a(E)/E = A © ^ © and the noise term (C) is 
consistent with zero for both the data and Monte Carlo 
and is removed from the fit. The sampling term B for 
electrons is 0.4 99 ± 0.008(\/GeV) from data while it is 
0.504 ± 0.006 (VGeV) in Monte Carlo, showing extremely 
good agreement. 

Finally, a very important parameter of the calorime- 
ter is the difference between electron and hadron re- 
sponses as a function of energy. Section |l^ describes 
how the hadron resolution and non-linearity depend crit- 
ically on this difference. In other words, the more simi- 
lar the electron and hadron responses are, the better the 
calorimeter resolution, and the more linear. The mea- 
sured electron/hadron difference must be corrected by 
1% to account for the fact that the electrons in the cali- 
bration beam started upstream, while those from neu- 
trino interactions (or those from hadron showers) are 
much more uniformly distributed throughout each steel 
plate. The ratio of the reconstructed energy of electrons 
compared to that of hadrons is 1.06 ±0.03 at 75 GeV, 
which corresponds to a small non-linearity, in agreement 
with what is seen in the calorimeter response as a func- 
tion of energy. 



15 Conclusions 

In this paper we outline the calibration technique and 
subsequent testing of the NuTeV calorimeter. The overall 
gain and time dependence of the calorimeter are tracked 
using muons from the neutrino beam, and PMT gains 
are determined to better than one percent for a given run 
period. Using several sets of linear ADC's with different 
gains, we are able to cover the large dynamic range re- 
quired to reconstruct minimum energy deposition at the 
percent level, as well as reconstruct energy deposition 
from 600 GeV neutrino induced hadron showers. 

Although the technique of calibrating the detector 
with muons from neutrino interactions may seem simple 
and straightforward, checking this technique requires a 
very detailed and well-designed calibration beam. By us- 
ing a low mass spectrometer with long lever arm, NuTeV 
is able to achieve event-by-event momentum resolution of 
better than 0.1%. In addition, as a consequence of care- 
ful calibration of the magnets in the spectrometer and 
measurement of the particle composition of the hadron 
beam, the absolute hadron energy scale of the calorime- 
ter is determined to better than 0.3%. 

The non-linearity in the hadron response of the 
calorimeter is measured and agrees with predictions 
based on the measured difference between the hadron 
and electron response at a particular energy. Finally, 
the muon and electron responses of the calorimeter are 
shown to agree with a GEANT-based Monte Carlo pre- 
diction, once the details of the calorimeter geometry are 
accurately included. The vital statistics of the NuTcV 
calorimeter are summarized in Table 0. In conclusion, 
NuTeV has accomplished its goal of calibrating the ab- 
solute energy scale of its calorimeter to the level dictated 
by the physics analyses that NuTeV is performing. The 
calibration beam data also yields a wealth of informa- 
tion about hadron, electron, and muon interactions in 
an iron-scintillator sampling calorimeter. These can be 
used to study designs of future calorimeters with simi- 
lar geometries, such as the MINOS detector to search for 
neutrino oscillations, as well as foiLspace-based calorime- 
ters to measure cosmic ray fluxes a. 
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Table 7: Vital statistics of the NuTeV calorimeter. 







Hadron Non-linearity 
from 5.9 GeV to 190 GeV 


3.0 ±0.5% 


naLLi oil ijjiit-ijiy ocdit; uiicciLaiiiLy 


u.^O /o 


Hadron Energy Resolution: 

/t(' P — i tn 
cr[-tj)/-tj — ^ 

A: Constant Term 

B: Stochastic Term {VCeV) 


0.022 ±0.001 
0.86 ±0.01 


Residual position dependence of 
hadron energy reconstruction 

... „ jl CA ... ^ ^^^^ 

more than 51) cm from edge 


^ n CO/ 
< 0.570 


Electron/Hadron Response Ratio 

(iminp' (Trrnnm PjirPimptpriy.fitin'n 1 


1 08 ± 01 


Electron Energy Resolution: 

A: Constant Term 

B: Stochastic Term {^/GeV) 


0.042 ±0.002 
0.499 ± 0.008 


Average Number of 
Photoclectrons /counter/MIP 


30 


Hadron MIP-to-GeV Conversion 
Factor (C^ GeV/MIP) 


.212 ± .001 



and Alignment groups. 
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